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Abstract 

Using the results of several quenched lattice simulations, we predict 
the value of the strange and charm quark masses in the continuum at 
the next-to-leading order, m^'^^(/i = 2 GeV) = (128 ± 18) MeV and 
m^(/i = 2 GeV) = (1.48±0.28) GeV. The errors quoted above have 
been estimated by taking into account the original statistical error of 
the lattice results and the uncertainties coming from the matching of 
the lattice to the continuum theory. A detailed presentation of the 
relevant formulae at the next-to-leading order and a discussion of the 
main sources of errors is also presented. 
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1 Introduction 



Quark masses are parameters of the QCD Lagrangian that cannot be deter- 
mined within QCD by theoretical considerations only and cannot be mea- 
sured directly since quarks do not appear as physical states. They are how- 
ever very important for several reasons. In the framework of GUTs, quark 
and lepton masses at low energies are related to the pattern of symmetry 
breaking at the grand unification scale. The best-known example is that, 
within GUTs, one predicts nib ~ 3m^ as a consequence of the SU{5) re- 
lation rrib = rrir at the GUT scale and of strong-interaction effects 
Constraints on quark masses and relations between quark masses and ma- 
trix elements of the Cabibbo-Kobayashi-Maskawa mixing matrix are also 
found by using some ansatz on the form of the quark mass matrix 
In this framework, for example, the matrix element \Vcb\ depends, among 
other things, on the ratio of the strange to bottom quark mass. The mass 
of the strange quark is also a crucial factor appearing in the evaluation of 
matrix elements of penguin and electro-penguin operators, which enter in the 
A/ = 1/2 K —>■ nir amplitude and in the CP violation parameter e'/e. The 
values of the light-quark masses are necessary also to estimate the quark- 
antiquark condensate (OlV'V'lO) which is the chiral order parameter of the 
QCD vacuum and finally the value of the up quark mass is fundamental for 
our understanding of the strong CP problem. 

Lattice QCD is in principle able to predict the mass of any quark by 
fixing to its experimental value the mass of a hadron containing a quark 
with the same flavour. The quark mass that is directly determined in lattice 
simulations is the "bare" lattice quark mass m{a), which can be converted 
to the continuum, renormalized mass m{fi) through a well-defined procedure 
[B^lUni- '^^^ conversion factor relating m(a) to m(/x) can be computed in 
perturbation theory. Since the typical scale is of order 1/a (or /i), where a 
is the lattice spacing, and in current numerical simulations ~ 2-4 GeV, 
we expect small non-perturbative effects. This would imply an accurate 
determination of quark masses at a scale larger than a~^. In this respect 
lattice QCD is unique, since other techniques, like QCD sum rules 



or 



have to work at much smaller scales, where higher-order corrections |T7 
non-perturbative effects [0 may be rather large. On the lattice, however, 
at values of the lattice strong coupling constant P = 6/g1{a) currently used 
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{P = 6.0 — 6.4), perturbative corrections are often quite sizeable, due to the 
presence of "tadpole" diagrams, which are absent in continuum perturbation 
theory, and may give rise to large uncertainties. In order to reduce these 
uncertainties, two methods, based on "boosted" perturbation theory [1^, ^ 



or on non-perturbative techniques |^ , have recently been developed and will 
be used in the following. 

In this paper, using values of quark masses determined in several nu- 



merical simulations we predict the corresponding masses in the 

continuum and discuss the uncertainties entailed by the use of perturbation 
theory at the next-to-leading order. Following ref. [^, we also envisage a 
procedure to determine the continuum quark mass, which completely avoids 
the use of lattice perturbation theory. Our predictions are essentially limited 
to the strange and charm quark masses. Actually, simulations do not have 
direct access to very light (up and down) quarks, because of finite-volume 
effects. Moreover one expects that the "quenched" approximation can intro- 
duce uncontrolled systematic errors in this case. Thus we do not have much 
to say about the up and down quark masses. On the other hand, we cannot 
put the b quark on the lattice, because it is necessary to satisfy the condi- 
tion ruba -C 1 in order to avoid huge discretization errors. However, mb can 



be obtained from quarkonia spectroscopy p% 33 1. An alternative method, 



based on the heavy-quark effective theory on the lattice will be presented 



elsewhere fB3 



In comparison with QCD sum rules, which is at present the only alterna- 
tive approach to predicting quark masses, lattice QCD has advantages and 
disadvantages. The main advantage is that, on the lattice, one is able to 
work at scales large enough to avoid higher orders or non-perturbative ef- 
fects which plague QCD sum rules calculations. Moreover, by making use of 
the proposals of refs.[]TU|, |2D| or |]2T[], the error due to the large perturbative 



corrections encountered on the lattice can be reduced. Finally it is possible 
to reduce discretization errors, by using some "improved" lattice quark ac- 
tion, following the proposal of Symanzik |^, |3^ . The main disadvantage 
of the lattice approach is that essentially all the results have been obtained 
in the "quenched" approximation and it is difficult to quantify the size of the 
error introduced by this approximation. 

We beheve that the results of this work are useful as a presentation of the 
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method, for the discussion of the uncertainties coming from the matching of 
the lattice to the continuum theory and for phenomenological apphcations. 
The plan of the paper is the following. In sec. ^ we give the formulae 
necessary to relate matrix elements of lattice operators and masses to the 
corresponding quantities in the continuum. The formulae introduced in sec. 
1^ are applied to quark masses in sec. ^ where lattice perturbation theory and 
non-perturbative methods based on Ward identities are reviewed. Numerical 
results and estimates of errors are given in sec. |[ 



Bare and renormalized operators on the 
lattice 



In this section we summarize the main formulae that relate lattice operators 
to the corresponding continuum renormalized operators^. These formulae 
will be used to give the continuum renormalized quark mass in terms of the 
lattice bare mass, as computed in Monte Carlo simulations. 

Let us consider the renormalization of a generic lattice two-quark oper- 
ator of the form Or (a) = V^r?/', where F is one of the Dirac matrices and 
ip is the bare lattice fermion field. At the next-to-leading order (NLO) the 
generic, forward, two-point Green function, computed between quark states 
of virtuality -C has the formQ 



r(pa) 



gMfl pa 
167r2 ^2 V ^ 



+ i(7"' + (-2,8„ + 7'°')C' + ,<?;(A^))ln(^)^)]+,,, (1) 



dX 



where . . . represent terms beyond the next-to-leading order and terms of 
0(a), which we will assume negligible in the following; Tq is the zeroth- 
order Green function, and g1{a) = 6/(3 is the bare lattice coupling con- 
stant; 7^°-* and 7^^-* are the leading (regularization- independent) and next- 



^ A more detailed discussion on this point can be found in ref. 
^ We work on a Euclidean lattice. 
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to-leading (regularization-dependent) order anomalous dimensions, respec- 
tively; A = A (a) is the lattice gauge parameter of the gluon propagator 

(llij,,y{q'^) = —6^u + (1 ~ ^)(ln(lv I (f^ ■ It obeys the renormalization group 
equation: 

1 ^?M = _,,(<,i(„)) = _,30?M+..., /3« = 5i^, (2) 

A(a) oa loTT^ ott 

We have introduced a scale-dependent gauge parameter in order to define 
a gauge-independent anomalous dimension and simplify the renormalization 
group equation for T{pa), see eq. (^) below. 

The lattice coupling constant obeys the equation: 

- -PiaM) - 2/5o^ + 2/?ip^, (3) 
where /3o,i are given by: 



{llN-2nf) ^ 34^. 10,, {N^ - I) 

and Uf is the number of flavours. Equation (|l|) guarantees that all the matrix 
elements can be made flnite (as a ^ 0) by multiplying the bare operator by 
a suitable renormalization constant, obtained by flxing the renormalization 
conditions for Or- 

From the above equations we flnd: 

la^ - Pialia))^^^ - Px{9l{a))X^ - i{gl{a)))T{pa) = (5) 
with 

^r„2 ( ^^ _ alia) (0) atia) (1) /p^ 

In view of the comparison with some continuum regularization, it is con- 
venient to expand the bare Green function T{pa) of eq. ([1|) in terms of the 
continuum minimal subtraction (MS) coupling constant, evaluated at the 
scale vr/a. The continuum MS coupling a<j(7r/a) is related to the lattice bare 
coupling (a) = g1{a)/4:7i by the equation: 

^ ^ ■ "^^"/^U + ...), (7) 



(a) asiTc/a) ^ An 
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where A is a numerical constant. With this substitution, eq. becomes: 



d 



d d 
ot.-^)- /5A(as)A— - 7L(as))r(pa) = 



das 



and 



(47r)^ 



a, 
47r 



(47r 



(8) 



(9) 



(10) 



By changing the expansion parameter, one also has to change the two loop 
anomalous dimension |^-|^, see also eq. (|T]): 



„(i) 



11' = T^' - 

Finally the running coupling constant is given by: 



4:71 



1 



QCD. 



(1 



QCD) 



+ ■ 



(11) 



(12) 



The above equation defines the continuum MS scale parameter Aqcd at the 
NLO. 

In continuum dimensional regularizations, after subtraction of the poles 
in 1/e, the renormalized Green function has a form similar to eq. (|I|): 



r(^) = To 



1 + 



a^(/i) 1 



i^(0)in(^)2 + c'^) 



+ 



It obeys the renormalization group equation: 



jj' 2 

P^2 



+ 



d 



d 



d 



(/i^ + /5(a,)^ + ^x{as)\^ + 7(«s))r(^) = 



da 



(13) 



(14) 
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with 



and 



Air 



^(0) + 



a,. 



(47r)2^ 



(1) 



1 dm n, ^ 

— = P\[as)- 



(15) 



(16) 



' A(yu) (9/i 

By solving the renormahzation group equations for T{pa) and r{p/fi), and 
imposing the matching condition V{p/ ji) = C{fia)T{pa), we find [0^] : 

Or(/i) = C(/ua)Or(a) 

«s(/i) x7(°V2/3or , a,(/i) -a,(7r/a)7W/?o-7^°^A 



+ 



Vas(7r/a)^ 



1 + 



47r 



47r 



Or(a) 



/ a,(/i) N7<°V2ft./ a,(/i) - a,(7r/a) 7(^)/jo - 7W/?i 



X 



«s(7r/a) 



1 + 



Or(a) 



(17) 



Equation ([T7| ) has been obtained using the universahty of the combination^: 
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(1) 



2Po 



C 



2Po 



(18) 



Equation (0) is interpreted as follows: the lattice operator is matched into 
the continuum operator at the scale vr/a through the factor 1 + agiii / a) / An 
— C^^ and then evolved in the continuum, according to eq. (0), from 
the scale vr/a up to /x. A few comments may be useful at this point: 

• One can eliminate the coefficient in eq. (|l7D by defining a new 
lattice operator 



Of (a) = Or(a)(l 



47r ' 



(19) 



This relation is true only when we use the same coupling constant and renormalized 
gauge parameter on the lattice and in the continuum. 
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This is equivalent to the regularization independent definition of the 
renormahzed operators discussed in refs.||38||-[40[]. This definition is 
such that r(p = 7r/a)/ro = 1 + O(a^)0. Notice however that this 
procedure requires a knowledge of the external states for which we 
have computed and that it is in general gauge- dependent. 

We have decided to expand the lattice Green function in ^^(Tr/a). How- 
ever we have the freedom to expand in ^^(l/a) or any other scale we like 
since the change will be completely compensated at the NLO. For the 
same reason we can expand the Green functions in a different coupling 
constant, for example the "boosted" coupling defined in refs.|]T^, EO]. 



There, it was argued that the series in may minimize 0{a1) NNLO 
corrections. If we expand in we have accordingly to reorganize the 
expression used for the coefficient function. It is wise, however, to check 
the stability of the physical amplitudes under a change in the scale used 
for as and assume this as a theoretical uncertainty. 



3 Quark masses 

3.1 Standard perturbative approach 

From lattice simulations, by fixing the value of the lattice spacing and the 
physical mass of a strange or charmed hadron, it is possible to evaluate the 
bare lattice quark mass m(a) = ms{a) or mch{cL) (strange and charmed quark 
masses). From m(a) one can compute the quark mass m^'^^lfi), renormal- 
ized in the continuum minimal-subtraction dimensional scheme, at the NLO; 
m^^{n) will be defined below. In the following we report the relevant for- 
mulae used to relate m^^{fi) to m(a) and discuss the different sources of 
theoretical uncertainty. The final values and errors of m'^^^{^) for different 
flavours have been obtained by using the results of refs.||22|-|3l|, as shown in 
table ^ and include the theoretical uncertainties discussed below. 



^ This is not in contrast with the condition ^ [tt / of' that we have to impose in order 
to avoid discretization errors. The renormahzation condition is simply a consequence of 
eq. (|l|), which is derived by expanding in as(7r/a) for <^ {71/0)^ . 



For definiteness, the formulae used in this section are valid for the lat- 



tice Wilson [g3] and SW-Clover ||35| quark actions. The case of staggered 
fermions will not be considered here. 

In perturbation theory, the inverse quark propagator on the lattice can 
be written as: 



S ^{p, mo) 1- (pa, 



moa) 



+ mo 



1 — S2 (pa, moa) 



+ 



yL 
^0 



(20) 



where a is the lattice spacing. The last term in eq. (pOf ) diverges linearly as 
a — i> 0. This term is induced by the explicit chiral symmetry breaking that is 
present in the lattice formulation of the quark action. The linear divergence 
is a mass term, which can be eliminated by a suitable redefinition of the bare 
quark mass 0-|lO[, 0- At 0{as), this is obtained by writing S'^{p, mo) 
as follows: 

S^^{p,m(a)) = 1 — Ti^(pa,m(a)a) 
X + m{a)(l + Jl^{pa,m{a)a) - E2{pa,m{a)a)^ , (21) 

where m{a) = m^ + Sg/a. By comparing eq. ( ^Ij) to the corresponding 
expression in the continuum: 



5-i(p,m^^(^)) 



Sf(p//i,m^^^(/i)//.) 



MS/ 



X 



ij^ + m*^^(/i)(l + Sf (p//i,m^^(/i)//i) - S2^(p//x,m^^(/i)//i) 



(22) 



we obtain: 



m^^^iyjj) = m{a) 1 + 'Z^(pa,m{a)a) — 'Z2{pa,m{a)a 

-Sf (p//., m^(/i)/^) + S^(p//., m^(/i)/^); 

~^(a) [1 + ^(7(0) ln{—) + K„ 
L 47r ^ jja 



(23) 



7^°'' and Km are numerical constants, which can be computed from the general 
expressions of Si 



Sf'^(p//i,m//i) 



as 



1 



471 2N 
a, iV^ - 1 
4^ 2N 



L,C 
^1 



0-2' + 



+ j rfx(2(l -x)ln((p^x(l -x) +m^x)//i^)) 
dx(4:\n{{p'^x{l - x) +m'^x)/fi^)'^] (24) 
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Action 


P 


a-' 


Kc 


Ks 


ms{a) 




Wilson PI 


6.0 


2.23(5)'' 


0.15704(3) 


0.1548(1) 


0.104(2) 


0.127(3) 


Wilson ^ 


1 


6.0 


2.21(8)" 


0.15704(3) 


0.1547(2) 


0.105(4) 


0.128(4) 


Wilson @ 


6.0 


2.29(37)" 


0.1570(1) 


0.1553(8) 


0.080(40) 


0.098(49) 


Clover [0 


6.0 


2.05(6)" 


0.14551(3) 


0.1437(1) 


0.090(3) 


0.103(3) 


Clover O 


6.0 


2.11(11)" 


0.14551(4) 


0.1438(2) 


0.088(4) 


0.100(5) 


Clover [g 


6.0 


1.95(7)" 


0.14548(2) 


0.1434(1) 


0.096(3) 


0.108(4) 


Clover [H 




6.0 


1.78(9)" 


0.14548(2) 


0.1430(2) 


0.104(5) 


0.117(6) 


Clover 






6.0 


2.01^ " 


0.14556(6) 


0.14371^ 


0.089(28) 


0.101(31) 


Wilson [|26 


1 


6.2 


2.9(2)" 


0.15330(2) 


0.1517(3) 


0.101(9) 


0.121(10) 


Wilson [|26 


1 


6.2 


3.0(2)" 


0.15330(2) 


0.1518(2) 


0.098(8) 


0.118(9) 


Wilson |27 


1 


6.2 


2.771^3" 


0.15328^4 


O.I517I3 


0.094(20) 


0.113(24) 


Clover [H 




6.2 


2.7(1)" 


0.14315(2) 


0.1419(1) 


0.082(8) 


0.093(8) 


Clover [p9|| 


6.2 


3.05(19)" 


0.14313(2) 


0.1421(1) 


0.079(5) 


0.090(5) 


Clover 


6.2 


2.73(17)" 


0.14313(2) 


0.1418(2) 


0.089(6) 


0.100(7) 


Wilson @ 


6.3 


3.01(26)" 


0.15163(3) 


0.1503(3) 


0.088(21) 


0.105(25) 


Wilson g 


6.4 


3.7(1)" 


0.1506(1) 


0.1495(1) 


0.089(3) 


0.106(3) 


Wilson 


6.4 


4.0(6)" 


0.1506(1) 


0.1497(2) 


0.082(13) 


0.098(15) 


Action 


[3 










rfichia) 


Wilson p6 


1 


6.0 


2.18(9)" 


0.15702(4) 


0.1313(16) 


1.36(4) 


1.26(2) 


Wilson p6 


1 


6.0 


1.97(8)" 


0.15702(4) 


0.1271(17) 


1.48(5) 


1.31(2) 


Clover 


M 




6.0 


1.92(11)" 


0.14544(2) 


0.1212(24) 


1.32(8) 


1.13(3) 


Clover [0 




6.0 


1.94(5)" 


0.14544(2) 


0.1218(11) 


1.30(4) 


1.12(2) 


Wilson ^ 


1 


6.2 


2.9(2)" 


0.15330(2) 


0.1354(24) 


1.24(8) 


1.22(5) 


Wilson |6 


1 


6.2 


3.0(2)" 


0.15330(2) 


0.1362(22) 


1.22(8) 


1.21(5) 


Wilson ^ 


1 


6.2 


2.771^3" 


0.1532811 


0.1350(-) 


1.22(10) 


1.20(10) 


Clover [|1 




6.2 


2.6(2)" 


0.1431311 


0.1290(-) 


0.99(8) 


0.94(7) 


Wilson |0 


1 


6.4 


3.7(1)" 


0.1506(1) 


0.1386(7) 


1.06(2) 


1.10(2) 


Wilson |0 


1 


6.4 


4.0(6)" 


0.1506(1) 


0.1402(24) 


1.00(11) 


1.05(9) 



Table 1: Lattice results for the strange and charm quark masses from 
lattice QCD. In the table (3 = 6.0/g1{a), Kg (q = s,ch) is the lattice 
hopping parameter for the quark q, amq{a) = l/2{l/ Kq — and 
afhq{a) = \n{4:Kc/Kg — 3). , m{a) and rh{a) are given in GeV. The 
symbol " denotes the scale taken from Mp and " from In the second case, 
we used boosted perturbation theory for the renormalization constant Za in 
the Wilson case (Za = 0.78,0.79,0.11 at /3 = 6.0, 6.2, 6.4^, and the non- 
perturbative values in the SW-Clover case (Za = 1.09,1.04 at (3 = 6.0,6.2 
Ref. ^ used Za = 0.77,0.79 at (3 = 6.0,6.3, respectively. 



Regularization 


0-1 


0-1 


0^2 


yL 


MS 


1 




2 




Lattice- Wilson 


13.85 


0.99 


1.90 


51.43 


Lattice-Clover 


10.11 


2.12 


-8.07 


31.98 



Table 2: Quantities entering in the definition of Km, eq. (]^^; Sq , defined in 
eq. ^^), is given in units of {as/4:7r){N^ — 1)/2A^. They have been computed 
m refs.^-^^,g^,^. 



L C r— 1 r 

where the constants 2 reported in table together with Sq , which 
was defined in eq. (pUj). One finds: 



7(0) = e: 



2N 



and 



iV^-l 



2A^ 



c 



(25) 

(26) 



1 1 — (T2 — 0"i + — 6 ln(7r) 

From eq. (^) we can derive the relation between m^^^{fi) and m(a) at the 
next-to-leading order, cf. eq. (|r^) in sec. 0: 



MS 



(/^) = ( 



as(7r/a) 



^ ^ a.(/i)-a.(7r/a) ^7W _ 7^°^^ ^ 



47r 



2/5o 
a^(7r/a) 



47r 



m(a); (27) 



■m[a), m^^dx), 7*^^^ and are regularization-dependent, gauge-invariant 
quantities; 7^°-* is the LO anomalous dimension and 7^^-' the NLO one, com- 
puted in the MS dimensional scheme; 7(1) is connected to the lattice anoma- 
lous dimension 7];^^ through the relation, see also eq. (p^Sj) : 



7(^)=7?+2/5oir„.. 



In ref. [^-|^ they found: 



7^ 



97^A^2 _ 

y 2N 



1 

- + 3 



lOn/ 



2N 



2N 



(28) 



(29) 
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In eq. ([27|), is the running coupling constant in the continuum minimal- 
subtraction scheme. It is related to the lattice bare coupling af'(a) = g1{a)/4:7i 
by the equation, see (J^): 



(a) <ys{q) ^ 47r qa 
At (3 = 6.0 one finds: 



as{q = ^) - 1-45 (1/a). (31) 

The above definition gives values of as close to the values of the "optimized 
perturbative lattice expansion" couplings aY introduced in refs.[|^, One 
typically finds a^(7r/a) ~ 1.6-1.8 (a). Among the definitions used in 
refs.flTSl, we have: 



1 1 ^ , o^'fe) 



where aj^i is defined from the Q-Q heavy-quark static potential. Using the 
perturbative expansion of the plaquette, one finds: 

_±_ = I (i + «r^(g = ^/Q) 6 ^33^ 

«f(a) «r^(g = 7r/a)(|Trf/,,,,)V ^ 47r 'J- ^ ^ 

One can expand in a^^ by redefining 7'-^'' according to eqs.(^ and (0), and 
use the value of a^^, found by using the plaquette computed in numerical 
simulations. Another commonly used definition of is: 

' -<«^^'=>'- + (34) 



V2 



We will discuss the uncertainties in the quark mass, coming from different 
choices of the expansion parameter in sec. ^ 

We can give an alternative definition of the renormalized mass, which is 
regularization- independent. One can for example define the renormalization 
constant of the mass as the inverse of the renormalization constant of the 
scalar density S = ipip- 

Zm = Zg , (35) 
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where Zs is defined by the renormahzation condition: 



Zs{iJ^a)rs{^a)/r 



(36) 



imposed to the amputated Green function on off-shell quark states with 
= yU^. Since Zs in eq. (|36|) is gauge-dependent, we will denote it as 
^gauge-ianjey,... ^-^^ foUowiug. The reuormalizatiou scheme introduced in 
eq. (^), which we call in the following RI for regularization-independent, is 
particularly suitable for the implementation of the relevant Ward identities 
of the theory, see sec. |3.3|, and for the non-perturbative renormahzation of 



lattice operators [^. In perturbation theory, in the Feynman or Landau 
gauge, one finds: 



fey, Ian 



\ (y A TT I n\ ^ 



+ 



'as(vr/a) 
asin/a) 



a^(/i) 



Ian 



r<fey 



Air 47r 
An ™ 



An 
m{a) 



2Pi 



ian 



2N 



2N 



-A 



(37) 



As explained in ref. |^0|5 ?Tj,-^^^''""(/i) are regularization-independent. In fact 
they are obtained by imposing the same renormahzation conditions in all the 
regularizations. The price is that this definition can be gauge-dependent and 
one has to specify the external states at which the renormahzation conditions 
have been imposed. 



3.2 Tadpole resummation 

In the Wilson or SW-Clover lattice formulation of the quark action, m(a) is 
expressed in terms of the hopping parameter as: 

m(a) =lnfl + -(— -^)) ~ -f— -^V (38) 
^ ' ^ 2^K K^J) 2VK K^J' ^ ' 
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where the critical value of K, corresponding to m(a) = 0, is K'^ = 1/8. At 
first order in perturbation theory, is given by: 

^ = 8 - = 8no. (39) 

It has been argued that lattice perturbation theory is dominated by tadpole 
diagrams, which have no correspondence in the continuum, and that the 



hopping parameter should be replaced by an "effective" one [|1^, |20 

K 



K = Ku> 







SK, 



(40) 



where Kc is defined non-perturbatively as the value of K at which the pseu- 
doscalar meson mass vanishes. The same tadpole diagrams contributing to 
Sq also enter the calculation of Si. At first order in perturbation theory, we 
can write: 



(41) 



1 - = 


yL 


-(Ef 


yL 
4 




2f) = 


8/C<^ 


-Sf) 



Using eq. (0), the lattice inverse quark propagator becomes: 

1 — Sf'(pa, m(a)a) 



S ip,m{a)) 



+ m{a) (l + (pa, m(a)a) — E2 (pa, m(a)a)^ (42) 

where m(a) = 4Kc/K — A ^ ln{AKc/ K — 3). The relation between the lattice 
quark mass m(a) and m*^'^(/i) is then given by the same expression as in eq. 
), with rh{a) instead of m(a) and cri = o"! — Sq/4 instead of ai (see table 



I)- 



3.3 Renormalization conditions and Ward identities 

In this section we discuss the definition of the quark mass via the Ward 
identities of the regularized theory and its relation with the quark mass 
defined in the previous section^. 

^ To avoid problems connected with anomalous terms, we will only consider Ward 
indentities for non-singlet currents. 
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In the continuum, the renormahzation of the quark mass and of the scalar 
density are related by the vector current Ward identity: 



m4 



m (43) 



where the ip^s are bare fields, M is the bare mass matrix and are 
arbitrary on-shell physical states. Equation ( ^31 ) guarantees that the product 
mipip is unrenormalized by strong interactions: 

mnijRi^R = mZmZsi^ip, (44) 

i.e. that Zm = Z^^ as in eq. (pSj), which remains valid on the lattice, where it 
is also possible to define a conserved vector current in the limit of degenerate 
quark masses. It ensures that the mass and scalar density renormahzation 
constants are the inverse of one another to all orders in perturbation theory 
and beyond. Still we have the freedom to decide which definition of the 
renormalized mass we wish to use: renormalized on the mass-shell, in the 
MS (gauge-invariant) or in the RI (gauge-dependent) schemes, as explained 
in sec. ^ Of course we could decide to renormalize the mass in the minimal- 
subtraction scheme and the scalar density in RI. Such an exotic choice 
would only obscure the understanding of the Ward identities and will not be 
considered here. Differences between different definitions of the quark mass 
are related by terms computable in perturbation theory. This is true for 
heavy quarks with the renormahzation on the mass shell and for heavy or 
light quarks with the MS or RI renormalizations, provided that fi and vr/a 
are much larger than Aqcd- 

The quark mass can also be defined through the Ward identity of the 
axial current: 

{a\d^A';\[3) = {a\d^^^,^,^m = j^, y | i,m- (45) 

In a regularization that preserves chiralityQ, as for example naive MS, the 
scalar and pseudoscalar (P" = ipX"- densities have the same renormal- 
ization constants, Zs = Zp, so that the quark mass defined from eq. (^) 

^ This regularization strictly speaking does not exist. 
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or p5| ) coincide. This is not true in continuum regularizations, which break 
chirahty, hke 't Hooft-Veltman or dimensional reduction in the continuum, 
or in the Wilson (SW-Clover) formulation of QCD on the lattice, where an 
explicit term that violates chirality is present in the action. We now explain 
how the quark mass can be defined through the axial Ward identity in the 
presence of an explicit breaking of chiral invariance. The argument is made 
with the lattice regularization, but is general. 

At zero order in as the axial current Ward identity on the lattice can be 



written as 10, 44 



{a\d^Al\(3) = 2m(a|Pl/?) + (46) 

in terms of bare fields and masses. We have simplified the Ward identity 
by taking degenerate quark masses; is the chiral rotation of the Wilson 
term, which is of 0(a) only in the free-field case (~ aipD'^ip). Because of 
interaction, develops terms of 0{g\/a) and g\ |jlO|, which renormalize the 
mass term, the pseudoscalar density and the axial current. Equation (^) 
becomes then: 

ZA{a\d^Al\P) = 2(m-m{m)){a\P'^\P) + (47) 

where Q is of 0(a), including strong-interaction effects. Close to the chiral 
limit, and neglecting terms of 0(a), eq. (|^) becomes: 

ZA{a\d^Al\P) = 2(m-m,)(l-^)(«|P'^|/3) 

= 2(m-m,)^{a\P''\(3), (48) 

where the chiral value of the mass is defined by the equation rric = rh{mc) 
and corresponds to K = Kc in the notation used in sec. ^ m — rric — 
1/2(1/ K — 1/Kc). The factor Zp/Zs appears since S and P' = Zp/ZsP 
(and not S and P [0) belong to the same chiral multiplet, because of the 
symmetry breaking present on the lattice. The same would be true with the 
't Hooft-Veltman regularization. 

Even though Z^ and Zp/Zs can be both computed in perturbation the- 
ory, an alternative, fully non-perturbative, definition of the quark mass can 
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be given. Using the Ward identities relative to the axial vector and scalar 
pseudoscalar densities we can determine both Za and Zp/Zs |^ 0, |5l], ^ 



in a non-perturbative way. Then we can measure the ratio 2p of the matrix 
elements of d^A'^^ and P*^ on a physical state; typically one uses: 

_ (0|gMgk(p = 0)) 
^ (0|P-|7r(p = 0)) ^ ' 

and then, by saturating the Ward identity in eq. (^81) we can find the lattice 
bare quark mass. This avoids ambiguities on the mass definition of the kind 
m(a) = A{Kc/ K — 1) ~ \n.{AKc/ K — 3), see sec. |3j^[]. In practice we define 
m = m{a) as: 

Zg 

m{a) = Zax — X p, (50) 

To find the continuum, renormalized quark mass, we can proceed at this 
point in two different ways. We can work in perturbation theory, and we 
essentially recover the results found in eq. (|27|) of sec. |3.2| . Alternatively we 
can renormalize non-perturbatively the pseudoscalar density by imposing, 
on quark states of momentum = fj,"^ and in a fixed (Landau) gauge, the 
renormalization conditions 



z's^^ifia)Tsifxa)/T^s = 4""(/"«)rp'(H/r?. = zl?"(Hrp(/ia)/r° = 1. 

(51) 

We have used the fact that the renormalization conditions imposed on S and 
P may be chosen in such a way as to satisfy the continuum Ward identities, 
i.e. Zp/Zs = Zp"'/Zg'"'. The continuum mass is then given by: 

m'""(/i) = Z^"(/ia)^(a) = x (Zj?"(/ia))"' x p, (52) 



We can then use eq. (|37D to obtain m'^^^{p) from m^""'(/i). The advantage 
of this procedure is that perturbation theory is used only to relate quark 
masses in the continuum, thus avoiding the large perturbative corrections 



present on the lattice pT| . The possible disadvantage is that p has to satisfy 



the condition Aqcd -C p <^ 1/a to avoid large higher-order corrections 

^ There remains another ambiguity of 0(a) in the definition of the time derivative in 
eq. (^). This ambiguity can be made of 0{asa) with improved actions. 
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Figure 1: Histograms of the strange and charm quark mass distributions in 
arbitrary units, computed in the MS scheme at the scale jj, = 2 GeV. They 
have been obtained by considering all possible uncertainties discussed in the 
text. In the case of the strange quark all the results of table H have been 
included. In the charm case, in order to reduce discretization errors, we have 
only considered the results obtained at [3 > 6.2. 

and discretization errorsQ. It turns out that, using the SW-Clover action 
at /? = 6.0, it is possible to choose a value of/i~l/a~2 GeV, at which 
good agreement can be found between the determination of Zp^/Z^"" on 
quark states and ZpjZg as computed using the Ward identities [p]] , 
with a relatively small discretization error. The systematic error due to 
discretization can be estimated to be of the order of ~ 10-15%. Since at 
present Zp" has been computed only at /3 = 6.0 and only with the SW- 
Clover action, the results of the non-perturbative method will be reported in 
the next section, but they will not be included in the final evaluation of m^. 

® One could imagine determining Z^^"^ non-perturbatively on the lattice, at very large 
values of /?, thus avoiding large higher-order effects. This however would not solve the 
problem because, in order to get a physical quantity from a matrix element computed in 
current simulations, one has in any case to evolve the operators down to a scale [i, smaller 
than the inverse lattice spacing used in the numerical calculation of the matrix element. 
Thus large higher-order effects would be present anyway. 
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Figure 2: Strange and charm MS masses and errors at jj, = 2 GeV from the 
results of table 

4 Numerical estimates of quark masses from 
the lattice 

4.1 Uncertainties of the lattice calculations 

We are now ready to give the results for the continuum quark masses. These 
results have been obtained by combining the values of the bare lattice masses 
given in table ^ with the conversion factor in eq. (0), computed with all 
its possible allowed variations, value of the scale, a^, etc. The theoretical 
predictions are subject to the uncertainties listed below: 

• Statistical error on m{a)a: this error is usually in the range 5-20% 
depending on the value of (3 = Q/g\{a) and on the statistical accuracy 
of the simulation. 

• Calibration of the lattice spacing: the uncertainty on the value of 

in physical units enters in two ways: 

i) when we convert the dimensionless quark mass m(a)a, which is di- 
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rectly accessible in numerical simulations, to m(a); 

ii) in the scale used to evaluate Q;s(7r/a) or a;s(l/a), etc. 

Only the first effect is important, since the second one gives a mild, 
logarithmic dependence in a. The value of a enters not only in the 
conversion of m{a)a to m(a), but also in the determination of Kg {Kchj 
from the mass of a strange (charmed) hadron. There are several meth- 
ods to evaluate the lattice spacing. The most popular ones are from 
the string tension, from the mass of the p, and from f^^. They usually 
differ from one another by ~ 10%. 

Choice of the expansion parameter: we can use a^, or ap, 
or any other coupling constant that one believes makes the perturba- 
tive series converge rapidly. The differences are of 0{a1), but may be 
important with the lattice regularization [0, We have evaluated 
the conversion factor with all the three possibilities listed above. In 
particular, with the MS choice of a^, we have evaluated either by 
using eq. (^) or by using the "unquenched" LEP results [^, i.e. 
we have computed as from eq. (|l^) with = (340 ± 120) MeV. 

Effects of 0{al): Values of the mass that differ by terms of order a1 
are obtained by using the two equivalent versions of eq. ([l^), but with 
different scales as arguments of a^. We have allowed the scale in to 
vary between 7r/a and 1/a, which seems to us sufficient to cover a large 
range of possibilities. 

Tadpole-improved definition of the mass: We have evaluated the 
continuum mass from m(a) defined as 1/2{1/K — 1/Kc) and from its 
tadpole improved expression ln{AKc/ K — 3) with the conversion factor 
changed accordingly. 

Discretization errors: In the final estimate of the mass we combine 
the results obtained with the Wilson and SW-Clover "improved" ac- 
tion, which suffer from different 0(a) effects. The comparison is useful 
because in the SW-Clover case, discretization errors are of ©(a^a) and 
have been found in some cases as small as 5%, to be compared to 30% 
in the Wilson case. 
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• Quenching: It is not possible to estimate and correct the error due 
to the quenched approximation. We have only applied the following 
procedure. From m(a) one can get m^^{TT/a), which is subsequently 
evolved to /i = 2 GeV. It is clear that m^'^^ {t[ / a) scales in a (and for 
consistency in /i) with a quenched anomalous dimension and /3-function. 
However, by imposing as{i^ / ay^^'^^^'^'^ = ^^(Tr/a)""'''"*^"''^'*'^'^, m^^^{7i/a) 
has been evolved also with the unquenched formulae and the difference 
has been taken into account in the final error. This procedure does not 
pretend to correct for the systematic error entailed by the quenched 
approximation, but it is consistent with the observation that most of 
the quenched and unquenched results look very similar after a suitable 
rescaling of the lattice coupling constant. 



4.2 Results for the strange and charmed quark masses 

By changing the value of the lattice spacing, the expansion parameter a^, the 
scale at which is evaluated, etc., and taking into account the statistical 
errors of the lattice simulations, we obtain a pseudo-Gaussian distribution 
of the value of the mass from which it is possible to evaluate the theoretical 
error, following the method explained in ref. see fig. 1. The continuum 
determination, at the scale fi = 2 GeV, of masses and relative errors from 
different lattice simulations are reported in figs. |l] and |^. One notices a nice 
agreement between different simulations performed at different values of the 
lattice spacing and with different lattice quark actions. No systematic trend 
in /i or dependence on the lattice action appears for the strange quark. In 
the case of the charm quark, a slight reduction of the value of the mass seems 
to appear at larger values of (3. This could be due to terms of 0{mchO')- In 
order to reduce the systematic error due to these corrections, in the final 
estimate of the value of the mass, we have taken only the determinations 
of rrich at /3 > 6.2. Combining all the results together, our best estimates, 
which is also reported in the abstract, are given by: 

mf^(/i = 2 GeV) = (128 ± 18) MeV 

mj^(/i = 2 GeV) = (1.48 ± 0.28) GeV. (53) 

For completeness we also try to estimate m^, with the SW-Clover action. 
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using the non-perturbative method envisaged in sec. p.3| . From a recent study 
at /3 = 6.0 they found |1: 



2p{K = Ks) = 0.071(5), p = 6.0. (54) 

From the non-perturbative lattice renormahzation of P"', ai f3 = 6.0 and 
(pa)^ ~ 1, they estimated pl| : 

Zj?"(/ia ~ 1) = 0.49(1) (55) 

from which, using eq. (p2|), one finds: 

m[^''{fi ~ 2GeV) = (150 ± 18) MeV, (56) 
or, from relation (|37D: 

mf^{l2 = 2GeV) = (141 ± 17) MeV, (57) 

in reasonable agreement with the results of eq. (|53|). In eqs. (|56| ) and (p^) 
only the statistical errors are reported. A further uncertainty of ~ 10 — 15%, 
due to discretization errors, must be taken into account. A more accurate 



analysis can be found in ref. |21 



5 Conclusion 

We have shown that lattice simulations of QCD can give accurate and reliable 
determinations of quark masses, using well-defined procedures. At present, 
it is possible to evaluate most of the uncertainties with the only exception of 
the systematic error introduced by the quenched approximation. The values 
of the strange and charm quark masses found in this work can readily be 
used for phenomenological applications. 
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